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Abstract 

The focus of this article is on recent results on ultra- 
high energy cosmic rays obtained with the Pierre 
Auger Observatory. The world's largest instru- 
ment of this type and its performance are described. 
The observations presented here include the energy 
spectrum, the primary particle composition, limits 
on the fluxes of photons and neutrinos and a dis- 
cussion of the anisotropic distribution of the arrival 
directions of the most energetic particles. Finally, 
plans for the construction of a Northern Auger Ob- 
servatory in Colorado, USA, are discussed. 
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1 Introduction 

The Earth's atmosphere is exposed to a flux of 
energetic particles from space. Their energies ex- 
tend from the MeV range to at least 10^" eV. The 
non-thermal energy spectrum follows a power law 
dN/dE oc E'^] the spectral index 7 has a value 
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around —3; it varies in certain energy regions in- 
dicating interesting possible changes in the com- 
position and in the acceleration and propagation 
processes. 

Diffusive shock acceleration is thought to be the 
basic mechanism that gives energy to charged cos- 
mic ray particles. The maximum attainable energy 
is proportional to the particles' charge and to the 
product of shock velocity, magnetic field and size of 
the acceleration region [5^ . Supernova remnants are 
candidates for the acceleration of galactic cosmic 
rays up to energies of several 10^^ eV. Of particu- 
lar interest is the study of cosmic particles at the 
highest energies observed so far. As cosmic rays 
of energy greater than 10^^ eV are not confined 
by typical galactic magnetic fields, it is natural to 
assume that those are produced by extra-galactic 
sources. The list of the very few viable candidate 
sources for such energies includes active galactic 
nuclei (AGN), radio lobes of PR II galaxies, and 
gamma-ray bursts (GRBs). However, the short en- 
ergy loss lengths indicate that particles of energy 
above 10^° eV should come from sources within a 
100 Mpc sphere. Astrophysical sources within our 
Galaxy are disfavoured. Por a recent review of as- 
trophysical sources, see [7j). 

More than 40 years ago Greisen, Zatsepin and 
Kuzmin realized that the interaction of protons 
with cosmic microwave background photons would 
result in significant energy losses. The energy spec- 
trum would show a fiux suppression above a thresh- 
old energy of about 6 x 10^^ eV, the GZK effect 
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Figure 1 : Energy loss lengths [5j of protons and nu- 
clei calculated for a redshift z = 0. Iron nuclei have 
larger or similar attenuation lengths to protons up 
to 3 X lO^o eV. 



[T]. Similarly, heavy nuclei are broken up due to 
photo-disintegration. A compilation of energy loss 
lengths of protons and heavier nuclei is shown in 
Figure[T]0. Since the energy loss processes for pro- 
tons and nuclei are different the evolution of com- 
position during propagation has to be taken into 
account. 

Alternative, non-acceleration scenarios for ultra- 
high energy cosmic rays have been proposed, such 
as the decays of topological defects or super-heavy 
dark matter. All of these models predict high fluxes 
of gamma rays and neutrinos. Finally there are 
propagation models in which the GZK energy loss 
processes are evaded or shifted to higher energies. 
Examples are violation of Lorentz invariance, the 
Z-huTst model or postulation of new particles with 
properties similar to protons. Reviews of such al- 
ternative scenarios can be found in [51 [S] . They are 
now disfavoured except for the highest energies due 
to already stringent limits on the flux of photons. 

Measurements of the arrival direction distribu- 
tion, primary mass composition and energy spectra 
are the keys to solving the puzzle of ultra-high en- 
ergy cosmic particles. High statistics is particularly 
important at the highest energies, where magnetic 
deflections are minimal and even charged particles 
are expected to point back to their source regions. 



2 The Pierre Auger Observa- 
tory 

The Pierre Auger Observatory for the highest- 
energy cosmic rays has been developed and built by 
more than 350 scientists in 17 countries. The south- 
ern site in Mendoza, Argentina, will be completed 
during the year 2008. The Observatory comprises 
1600 water-Cherenkov detectors deployed over 3000 
km^ on a triangular grid with 1500 m spacing. This 
surface detector array (SD) is overlooked by 24 elec- 
tronic telescopes, which are arranged in four sta- 
tions around the SD area. The telescopes record 
images of the faint, ultra-violet fluorescence light 
excited by the air showers. The fluorescence detec- 
tor (FD) operates on clear dark nights and achieves 
a duty cycle of 13 percent. The properties and per- 
formance of the prototype instruments have been 
published in 2004 [3] . A recent description of the 
combined hybrid operation of the SD and FD sys- 
tems of the Auger Observatory can be found in the 
proceedings of the ICRC 2007 [1|. 

The layout and deployment status as of June 
2008 is shown in Figure [21 The data set under- 
lying the results reported here has been collected 
basically from January 2004 to August 2007. Dur- 
ing this time the SD system grew from 154 to 
1200 tanks and the FD system from 6 to 24 tele- 
scopes. The exposure over this period is three times 
greater than that of AGASA, similar to the monoc- 
ular HiRes exposure, and corresponds to about 80 
percent of a full Auger-year. Above 1 EeV the 
Auger Observatory has recorded more events than 
all previous efforts together. Due restrictions to 
periods of good data quality have been applied. 
We propose to use the unit Linsley as a measure 
of the total exposure of cosmic ray observatories, 
ILinsley — Ikm^sryr. At the time of writing 
(June 2008) the total Auger exposure is 10,000L. 
The expected future annual increment is 7000 L. 

From the beginning the Auger Observatory had 
been conceived as a hybrid system, in which the 
SD and FD components are completely integrated. 
The merits of the SD include stable operation with 
100 percent duty cycle and a relatively straight- 
forward determination of the effective area times 
solid angle (the aperture), whereas the energy mea- 
surement has to rely heavily on Monte Carlo sim- 
ulations. Complementary, the FD offers optical 
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Figure 2: The southern Pierre Auger Observatory 
is located in the El Nihuel region in the Province of 
Mendoza, Argentina; the central campus is in the 
city of Malargue. Dots indicate the 1600 surface 
detector positions. The four named fluorescence 
detector stations with six individual telescopes each 
are all operational, as well as the surface detectors 
within the shaded area. 



shower detection in a calorimetric way and can be 
calibrated with very little dependance on shower 
models; however, the aperture grows with shower 
energy and its determination for all operating con- 
ditions is nontrivial. 

A public internet display of 1% of the events 
recorded by the Auger Observatory is available on- 
line at http://augerswl.physics.utah.edu/ED/, 

Each SD station is a 3.6 m diameter polyethy- 
lene tank containing a sealed liner with a reflec- 
tive inner surface. The liner contains 12,000 1 
of pure water. Cherenkov light produced by the 
passage of particles through the water is collected 
by three nine-inch-diameter photomultiplier tubes, 
which are symmetrically distributed at a distance 
of 1.20 m from the center of the tank and look 
downwards through windows of clear polyethylene 
into the water. The surface detector station is au- 
tonomously operating with a 10 Watt solar power 
system. The electronics package includes a pro- 
cessor, GPS receiver, radio transceiver and power 
controller. 

A vertically arriving cosmic ray shower of 10^^ eV 
typically triggers 8 detectors. SD signals are mea- 
sured in units of a relativistic muon passing verti- 



cally through the centre of a water-Cherenkov de- 
tector (vertical equivalent muon, VEM). This cali- 
bration of the SD array is performed continuously 
during data-taking with 2% accuracy. The number 
of secondary particles at observation level increases 
linearly with primary energy with some weak de- 
pendencies on the primary mass. The measured 
signal as a function of distance from the shower axis 
is fitted for individual events to obtain the VEM 
particle density at 1000m, S(IOOO), which is used as 
an estimator of the size of each event. For energies 
above 10^^ eV, the relative uncertainty of S(IOOO) is 
about 10%, including contributions from counting 
statistics, from the imperfectly known lateral parti- 
cle density distribution and from shower-to-shower 
fluctuations. 

At least five active stations must surround the 
detector with the highest signal, and the recon- 
structed shower core must lie inside an active tri- 
angle of detectors. For most analyses we use show- 
ers up to zenith angles < 60 degrees. These 
quality cuts result in full acceptance above 3 EeV, 
and nearly constant exposure as function of siri^O. 
The effective area of the entire array at any time 
has been calculated from the number of active 
hexagons, which can be deduced from the low- level 
triggers sent by each detector every second. The in- 
tegrated exposure is 7000 km? sr yr and is known to 
3%. Periods corresponding to less than 10% of the 
integrated exposure had data acquisition problems 
and were removed. 

The angular resolution of the surface detector 
was determined experimentally, checked using the 
pairs data set and found to be better than 2° for 
3-fold events {E < -i EeV), better than 1.2° for 4- 
fold and 5-fold events {Z < E <IQ EeV) and better 
than 0.9° for higher multiplicity events, which have 
more than 10 EeV. 

The FD measures the longitudinal development 
of cosmic ray showers in the atmosphere. Each in- 
dividual telescope of the FD images a portion of the 
sky of 30 degrees in azimuth and elevation. Light is 
collected by a segmented spherical mirror of 3.6 x 
3.6 m? through a UV-transparent filter window and 
a ring corrector lens to reduce the abberations in- 
herent in the Schmidt optics. The camera consists 
of 440 hexagonal photomultiplicrs, each with a field 
of view of 1.5 in diameter. The signals are contin- 
uously digitized at 10 Ms/s, temporarily buffered 
and searched for shower track patterns in real time. 



The timing of FD and SD systems are synchronized 
to about 120 ns. The data are merged offline in the 
event building process. 

The fluorescence light produced by the shower is 
detected as a line of pixels in the FD camera. The 
plane defined by the orientation of pixels, together 
with the timing information, is used to determine 
the shower direction. The pixel signal must be cor- 
rected for attenuation of the fluorescence light due 
to Rayleigh and aerosol scattering along its path to- 
ward the telescope, and for a contribution of direct 
and scattered Cherenkov light. The fluorescence 
light is then proportional to the electromagnetic en- 
ergy deposited by the shower along its path in the 
atmosphere. The reconstructed energy profile is 
fitted with a Gaisser-Hillas function [2] , which pro- 
vides a measurement of the maximum depth of the 
shower, and of the shower energy. A final correc- 
tion for missing energy due to high-energy muons 
and neutrinos amounts to about (10 ± 4)%. This 
procedure provides a nearly calorimetric, model- 
independent energy measurement with a resolution 
of 8%. 

The absolute fluorescence yield of the 337 nm 
band in air is 5.05 photons/MeV of energy deposit 
at 293 K and 1013 hPa, derived from [10]. We use 
the measurements of [11] for the wavelength and 
pressure dependence of the fluorescence spectrum. 

The absolute calibration of the FD telescopes is 
known to 10%. The relative response of all FD 
channels is determined twice every night illumi- 
nating the cameras from pulsed LEDs and/or Xe 
flashers. The absolute calibration is performed less 
frequently. In this case the whole aperture is illu- 
minated by a flat-field source with known spectral 
and directional characteristics and known intensity, 
calibrated at the National Institute of Standards 
and Technology. Atmospheric monitoring is an in- 
tegrated part of the FD operation. The system 
consists of four single wave-length LIDAR stations, 
one Raman LIDAR, several weather stations, bal- 
loon born meteorological probes, remote-controlled 
laser and dedicated aerosol monitors. 

The systematic uncertainties in setting the FD 
energy scale sum to 22%. The largest uncer- 
tainties are due to the absolute fluorescence yield 
(14%), the absolute calibration of the FD (10%) 
and the reconstruction method (10%), while the 
influence of pressure, humidity and temperature, 
the wavelength dependent response of the FD, the 



aerosol phase function, missing energy and others 
are smaller. 

Instrumental enhancements are currently being 
installed close to the Coihueco FD station. These 
include underground muon detectors, additional 
water Cherenkov detectors, high-elevation fluores- 
cence telescopes for a larger field-of-view and radio 
antenna to detect the geo-synchrotron emission of 
air showers. 



3 The Energy Spectrum 

Combining the SD and FD information, the en- 
ergy can be estimated for each event by a method 
that relies almost entirely on data and depends 
very little on hadronic interaction models or as- 
sumptions about the nature of the primary par- 
ticle. We use the so-called constant intensity cut 
method T2] , by which the zenith angle dependence 
of S(IOOO) is derived from the assumption that the 
true cosmic ray intensity at a given energy should 
be the same for all directions. Using a suitably 
chosen reference intensity we normalize S(IOOO) for 
each event to the signal it would have had at a 
zenith angle of 38° , which is the median zenith an- 
gle of the events of interest. That quantity called 
S'38 is used to correlate with the calorimetric en- 
ergy measurement by the FD for currently 661 se- 
lected high-quality hybrid events. To avoid effi- 
ciency and selection biases, low energy events are 
discarded from the calibration sample. The fit- 
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ted power law EpD = 1.49 x lO^^eF x S. 
dicates that S'38 grows approximately linear with 
energy. The energy resolution estimated from the 
root mean square deviation of the distribution is 
19%, which is in good agreement with the quadratic 
sum of the SD and FD energy uncertainties. 

The energy spectrum based on 20,000 SD events 
with the energy scale set by the FD as described 
above is displayed in Figure Sj The total system- 
atic energy scale uncertainty is 22% as described 
in the previous section. The residuals relative to a 
spectrum with a spectral index of 2.69 is also shown 
together with data from the HiRes experiment ^25j . 
No data from the AGASA experiment are shown as 
they are currently under revision |24j . The Auger 
data show the flux suppression above 4 x 10^^ eV 
with 6 standard deviations. 

Statistics and energy range can be extended by 
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Figure 3: Correlation between surface detector sig- 
nal and FD energy for 661 hybrid events used in 
the fit. The full line is the best fit to the data; 
the fractional differences between the two energy 
estimators are shown in the inset. 



including inclined showers with zenith angles larger 
than 60° and hybrid events, which have at least 
one SD detector together with a FD track recorded. 
The inclined events add another 1,500 L to the data 
set at all energies; special reconstruction efforts are 
currently being devised for events with zenith an- 
gles larger than 80°, where the particle distribu- 
tions are significantly distorted due to absorption 
in air and the Earth's magnetic field. 



4 The Primary Mass Compo- 
sition 

Measuring the composition of cosmic rays is crucial 
in order to obtain a full understanding of their ac- 
celeration processes, propagation and relation with 
galactic particles. The atmospheric depth Xmax 
denotes the longitudinal position of the shower 
maximum, which is directly accessible with the FD. 
It grows logarithmically with the energy of the pri- 
mary particle. The behaviour of Xmax for different 
primary particles like photons, protons and heav- 
ier nuclei can be conceptually understood in the 
framework of Heitler models [13' , which are in good 



Figure 4: Upper panel: Differential ffux as a func- 
tion of energy. Vertical error bars represent the 
statistical uncertainty only. The number of events 
in each bin is also given. Lower panel: Fractional 
differences between Auger and HiRes I data relative 
to a spectrum with an index of 2.69. 



agreement with detailed Monte Carlo simulations. 
Data and models for Auger hybrid events are shown 
in Figure [5l The appearent change from lighter to 
heavier particles above 2 x \0^^ eV is an interesting 
feature, which is currently under investigation. 



5 The Photon and Neutrino 
hmits 

Photon induced air showers would penetrate 200 - 
300 g/cvi^ deeper into the atmosphere than pro- 
tons. The observables sensitive to the longitudinal 
shower development are Xmax measured with the 
FD, and the signal risetime and the curvature of the 
shower front measured in SD-only events. The first 
limit on the photon contents of cosmic rays at the 
highest energies had been derived from the fluores- 
cence detector [T3] ; a more stringent result has now 
been derived using the large statistics of surface de- 
tector events. The maximum fraction of photons is 
2.0% above lO^^ eV and 62% above 8 x lO^^ eV, 
which corresponds to limits on the fiux of photons 



o 

n 

E 

X 

V 



- 










- 


^ Auger ICRC07 








'- 


p<°^°lssr;;^ 


-->" 


.--■¥" 


-"■■■■■'■■"'"t. 


•^ 




<^ 


1^'- 












QGSJETII-03 


- 


J7a ^,,--r:^^^-^---'"' 






----QGSJET01 


- 


^„.,''<^^^'.---'' 






SIBYLL2.1 


< 


-<^.-"' \<of\ 




, , 1 


-■-■ EPOS1.6 




10'^ 




10^^ 





E[eV] 



Figure 5: Atmospheric depth X^ax of the shower 
maximum as a function of energy. The data points 
are accompanied by their respective number of 
events in each bin; the expectations for primary 
protons and iron nuclei are shown for three differ- 
ent models 1^. 



of 6.9 X 10"'^ and 1.7 x Vd^^krrT^ sr~^yr~^ , respec- 
tively, with a confidence level of 95%. This new 
result 15J disfavors many exotic models of sources 
of cosmic rays as shown in Figure [6] . 

A large air shower array like Auger represents 
a significant target mass to high-energy astrophys- 
ical neutrinos. Neutrino induced showers can be 
identified if they occur deep in the atmosphere un- 
der large zenith angles, or by their special topol- 
ogy in the case of Earth-skimming tau neutrinos. 
Identification criteria have been developed to find 
EAS that are generated by tau neutrinos emerging 
from the Earth. No candidates have been found in 
the data collected between 1 January 2004 and 31 
August 2007. We derive an upper limit on the dif- 
fuse tau neutrino fiux as EldNy^/dEy < 1.3 x 10^'' 
GeV cm^^ s~^ sr~^ in the energy range 2 x 10^^ eV 
< E,j < 2x 10^^ eV. In Figure[7]we show our result 
|16j . which is at present the most sensitive bound 
on neutrinos in the EeV energy range. In the fu- 
ture, our sensitivity will improve by more than an 
order of magnitude. 



6 The Arrival Directions 

The search for anisotropics in the arrival direc- 
tions of cosmic rays or even for distinct sources 
has been a long-standing goal. The Auger Col- 
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Figure 6: The upper limits on the fraction of pho- 
tons in the integral cosmic-ray flux derived from 
SD events (black arrows) along with previous ex- 
perimental limits (HP: Haverah Park; Al, A2: 
AGASA; AY: AGASA- Yakutsk; Y: Yakutsk; FD: 
Auger hybrid limit). Also shown are predictions 
from top-down models and for the GZK photon 
fraction. See [H] and references therein. 
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Figure 7: Limits at 90% C.L. for a diffuse flux of 
astrophysical neutrinos. Limits from other experi- 
ments are converted to a single flavour assuming a 1 
: 1 : 1 ratio of the 3 neutrino flavours and scaled to 
90% C.L. where needed. The shaded curve shows a 
typical range of expected fluxes of GZK neutrinos 
- predictions vary by an order of magnitude. See 
[TO] and references therein. 



laboration has recently published the observation 
of a correlation between the arrival directions of 
the most energetic cosmic rays and Active Galactic 
Nuclei (AGN) listed in the Veron-Cetty and Veron 
catalogue [T7( ITS] . An exploratory scan of an early 
data set had been performed in the parameter space 
given by the angular deviation between arrival di- 
rections and AGN positions, AGN distance (given 
by their redshift) and cosmic ray energies. Param- 
eters were then fixed a priori and applied to an 
independent second data sample in order to avoid 
the complications and statistical penalties inher- 
ently present in a posteriori searches. The cor- 
relation has its maximum significance for energies 
greater than 57 EeV, AGN closer than about 71 
Mpc distance and angular deviations within 3.1°. 
In this conference contribution a different represen- 
tation in equatorial coordinates instead of galactic 
cooordinates is shown in Figure [5] 



This result has to be considered together with 
the flux suppression described previously; it is con- 
sistent with the hypothesis that the rapid decrease 
of flux measured by the Pierre Auger Observatory 
above 60 EeV is due to the GZK effect and that 
most of the cosmic rays reaching Earth in that en- 
ergy range are protons from nearby extragalactic 
sources, either AGN or other objects with a simi- 
lar spatial distribution. The maximum acceleration 
power of cosmic ray sources cannot be determined 
yet, but this alone cannot explain the sudden onset 
of the position correlation. 



A rich spectrum of interesting science questions 
opens up when future observations of all sources 
on the sky with high statistics can be combined 
with detailed investigations of the primary parcti- 
cles' nature and their interactions at center-of-mass 
energies up to 350 TeV. 



The Auger Collaboration has also searched for 
signals from the galactic center [TSj, and for clus- 
tering on different angular scales at the highest en- 
ergies and for correlations with BL Lac objects. 
These studies [501 IHl HI] have not confirmed pre- 
vious claims. 
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Figure 8: Map of the celestial sphere in equatorial 
coordinates. Colored circles with radius 3.1 degrees 
are centered at the arrival directions of the 27 cos- 
mic rays with highest energy detected by the Pierre 
Auger Observatory. The positions of the 472 AGN 
with redshift z less than 0.018 (corresponding to 
distances up to 75 Mpc) from the 12th edition of the 
catalog of quasars and active nuclei by Veron-Cetty 
and Veron are indicated by the blue stars. Grey 
shading indicates larger relative exposure, which is 
maximum at the South celestial pole. Centaurus 
A, Fornax A, and Virgo A, three of the closest can- 
didate sources, are marked with black stars. The 
super-galactic plane is shown as a dashed blue line. 

7 The Northern Pierre Auger 
Observatory 

The scientific challenge is to follow-up the recent re- 
sults by the identification of sources over the whole 
sky and by the detailed study of the physical pro- 
cesses at work in such extreme conditions. These 
involve strong gravitational fields, high tempera- 
tures and densities, large magnetic and background 
radiation fields, and very high energies interactions. 
Multiple events from the most intense sources will 
allow a detailed study of the source characteristics 
and acceleration processes, while particle collisions 
at center of mass energies in excess of 300 TeV 
may reveal new insights to particle physics. Such 
measurements at Auger South have yielded unex- 
pected results when compared with current interac- 
tion models. High statistics at the highest energies 
is essential to address these scientific objectives ef- 
fectively. The Northern Auger Observatory (Auger 
North) will focus on achieving higher statistics at 
energies above 6 x 10^^ eV, where the Greisen- 



Zatsepin-Kuzmin effect makes it possible to study 
nearby sources without the isotropic background 
from the rest of the far universe. 



Auger North will be built with a combination of 
the same basic elements as used in Auger South: 
surface particle detector stations, fluorescence tele- 
scopes, and an associated communications and cal- 
ibration infrastructure. The site, chosen in 2005, 
is located in the South-East corner of the State 
of Colorado (USA), centered at about 38° N Lat, 
102° 30' W Long. The average altitude is about 
1300 m above sea level. The landscape is almost 
flat, gently rolling and open; it offers an excep- 
tionally large area of more than 8,000 square miles 
(20,000 square kilometers), possibly spanning into 
the State of Kansas. An important feature is the 
system of county roads, which are spaced on a rect- 
angular one-mile grid covering a large fraction of 
the anticipated deployment area. 

The conceptual design includes a base grid of 
4,000 SD stations to be positioned on every sec- 
ond corner of the square-mile grid, alternating in 
adjacent lines. The total area is 8,000 miles^ or 
20,000 km^ with a detector spacing of approxi- 
mately \/2 X 1 mile (2.3 km). Simulations sug- 
gest that this conflguration reaches 50% efficiency 
at 10^^ eV and 100% efficiency at lO^^-^ eV. An 
embedded denser grid of 400 additional SD sta- 
tions would be placed on the empty positions in 10 
percent of the base grid in order to achieve full effi- 
ciency already at 10^^ eV for a subset of the events. 
The landscape requires a peer-to-peer data commu- 
nication rather than a direct transmission from SD 
stations to centralized receivers. Almost full cover- 
age of the SD system can be achieved with 40 to 
50 fluorescence telescopes anticipating a maximum 
viewing distance of 40 kilometers. 



Research and development works in order to im- 
prove performance and reduce cost have started 
some time ago and activities on the Northern site 
in South-East Colorado are ramping up. Techni- 
cally, the construction of Auger North could begin 
in 2010. Figure [9] shows the exposures of current 
and future cosmic ray observatories as a function 
of time. 
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Figure 9: Exposures of cosmic ray observatories 
as a function of time. The JEM-EUSO project is 
a space-born instrument on the ISS, whereas the 
other installations are ground-based. [26]. 



8 Conclusion and outlook 

The Pierre Auger Observatory has been collect- 
ing data of unprecedented quality since January 
2004. First results include a measurement of the 
energy spectrum, which exhibits a flux suppression 
above 40 EeV as predicted by Greisen, Zatsepin 
and Kuzmin. The arrival directions of the most en- 
ergetic cosmic particles are anisotropic and show a 
correlation with the positions of nearby extragalac- 
tic objects. Other previous claims of anisotropics 
or distinct sources have not been confirmed. Pre- 
liminary indications of the cosmic ray mass compo- 
sition give rise to interesting interpretations. Based 
on the wealth of current data, a conceptual design 
for the Northern Auger Observatory has been pre- 
sented. 
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